Combining directional light output and ultralow loss in deformed microdisks 
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A drawback of high-quality modes in optical microdisks is their isotropic light emission characteristics. Here 
we report a novel, robust, and general mechanism that results in highly directional light emission from those 
modes. This surprising tinding is explained by a combination of wave phenomena (wave localization along 
unstable periodic ray trajectories) and chaotic ray dynamics in open systems (escape along unstable manifolds). 
The emission properties originating in the chaotic ray dynamics permit directional light output even from mi- 
crolasers operating in the common multi-mode regime. We demonstrate our novel mechanism for the limafon 
billiard family and find directional emission with narrow angular divergence for a significant range of geometries 
and material parameters. 

PACS numbers: 42.55.Sa, 42.25.-p, 42.60.Da, 05.45.Mt 
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The possibility to confine photons in all three spatial di- 
mensions using microcavities has triggered intense basic and 
applied research in physics over thepast decade [JJ, e.g., re- 
search on ultralow threshold lasing QUI], single-photon emit- 
ters f7], solid-state cavity quantum electrodynamics fs',^,!?!]. 
Prominent examples of optical microcavities are whispering- 
gallery cavities such as microdisks flfefl, microspheres ifiol 
and microtori UH [ij] which trap photons for a long 
time T near the boundary by total internal reflection. The cor- 
responding whispering-gallery modes have very high quality 
factors Q = cjr, where uj is the resonance frequency. For mi- 
crodisks, the record Q -factor is around 6.6 x 10^ 1 9]. The high 
Q-factors and the in-plane light emission make microdisks at- 
tractive candidates for several optoelectronic devices, espe- 
cially for the nitride material system ifTsh where other cav- 
ity designs such as vertical-cavity surface-emitting laser (VC- 
SEL) micropillars faces severe challenges in mirror fabrica- 
tions lfr6ll . Unfortunately, the possible use of microdisks is 
limited by the fact that the in-plane light emission is isotropic. 

Shortly after the first fabrication of semiconductor mi- 
crodisk resonators it has been demonstrated that deforming 
the boundary of a disk allows for improved directionality 
of emission and therefore for more efficient extraction and 
collection of fight llZl [Ij, IjJ, |2Q]. Several shapes have 
been proposed and realized since then, but only few lead to 
light emission into a single direction with reasonable angu- 
lar divergence Jill that is essential for applications like 
single-photon sources. Moreover, all deformed microdisks 
discussed in the literature have a serious problem, Q spoil- 
ing [23]: The Q-factor degrades dramatically upon deforma- 
tion; in the worst case ruling out any application. The trade- 
off between Q-factor and directionality is not only a problem 
of microdisks but also of microspheres, microtori, and even of 
VCSEL-micropillars: for given total number of Bragg mirror 
pairs such a cavity can be optimized either w.r.t. directional- 
ity (significantly more mirror pairs at the bottom than at the 
top) or Q -factor (roughly equal number of mirror pairs at the 
bottom and top). 

Recently, a scheme to achieve highly directional emission 



without Q-spoiling has been suggested |24]. It employs the 
characteristic modifications of spatial mode structures near 
avoided resonance crossings. The drawback of this approach 
is the existence of nearly degenerate modes which can have 
different far-field patterns (FFPs). If two such modes are 
excited simultaneously (typical for present-day devices), the 
directionality might be lost. Nearly degenerate modes with 
similar FFP do occur but to find them requires laborious nu- 
merical calculations and a sophisticated adjustment of geome- 
try parameters depending on refractive index, wavelength and 
cavity size. 

In this Letter we introduce a novel, robust, and generally ap- 
plicable mechanism to achieve highly directional light emis- 
sion from high-Q modes in microdisks which does not suffer 
from the above-mentioned drawback. The key to our finding 
is to apply what is at the heart of quantum chaos |2?] and 
nonlinear dynamics flSn of open systems, respectively: wave 
localization along unstable periodic ray trajectories in systems 
with chaotic ray dynamics [i27i1 and the only recently acknowl- 
edged importance of the so-called unstable manifold for the 

FFPs of microcavities (ii m us m. The wave localiza- 
tion ensures the desired high Q-factors, whereas the unstable 
manifold provides the directional emission. Concerning the 
realisation of this scheme in practise, adequate microcavity 
devices can be expected to be easy to fabricate as no sophis- 
ticated adjustment of parameters is required. They are, more- 
over, well suited for multi-mode laser operation as all high-Q 
modes of given polarization possess similar FFPs. 

A microdisk is a quasi-two-dimensional geometry de- 
scribed by an effective index of refraction n. We assume 
n = 3.3 (GaAs) both for transverse magnetic (TM) and trans- 
verse electric (TE) polarization. A slight polarization depen- 
dence of n is neglected; it could be adjusted in the fabrication 
process, e.g., by changing the slab thickness. For the bound- 
ary curve of the deformed microdisk we choose the limagon 
of Pascal which reads in polar coordinates (p, ^) 



p((/)) = R{1 + ecosc 



(1) 



The corresponding family of closed cavities is known as 
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"limagon billiards" Ray and wave dynamics in bil- 

liards has been extensively discussed in the fields of non- 
linear dynamics ll26l] and quantum chaos ll25l] . The limit- 
ing case of vanishing deformation parameter e is the circle 
with radius R. Figure [Tta) illustrates a whispering-gallery 
ray trajectory in a circular microdisk trapped by total inter- 
nal reflection. A two-dimensional phase space representation, 
the so-called Poincare surface of section (SOS), is shown in 
Fig-ffib). Whenever the trajectory hits the cavity's boundary, 
its position s (arclength coordinate along the circumference) 
and tangential momentum sin x (the angle of incidence x is 
measured from the surface normal) is recorded. For £ = 0, 
rotational invariance of the system implies conservation of the 
angular momentum ex sin x- Ignoring wave effects, such a ray 
never leaves the cavity since it cannot enter the leaky region 
between the two critical lines for total internal reflection given 
by sinxc = 

Figures [Jc) and (d) show a trajectory in the lima9on cav- 
ity for £ = 0.43. In contrast to the case of small deforma- 
tion parameter s 1I23I1 the dynamics is predominantly chaotic. 
Starting with an initial x well above the critical line, a test ray 
(square, thick dots, and triangle) rapidly approaches the leaky 
region (sin x is not conserved) where it escapes according to 
Snell's and Fresnel's laws. Without refractive escape (n = 00, 
hard wall or closed billiard limit), the trajectory would fill the 
phase space in a random fashion (small dots). Periodic ray tra- 
jectories do exist but are always unstable, except for the two 
islands in the leaky region. Whispering-gallery trajectories 
are confined to the tiny region | sin x| > 0.99. 

While in the long-time limit the phase space of closed 
chaotic systems is essentially structureless, cf. Fig. [Hd), the 
phase space of an open chaotic system is structured by the so- 
called "chaotic repeller" [26]. It is the set of points in phase 
space that never visits the leaky region both in forward and 
backward time evolution. The stable (unstable) manifold of 
a chaotic repeller is the set of points that converges to the 
repeller in forward (backward) time evolution. The unsta- 
ble manifold therefore describes the route of escape from the 
chaotic system. In the case of light, Fresnel's laws impose an 
additional, polarization dependent weighting factor to the un- 
stable manifold in the leaky region |31], since at each reflec- 
tion the intensity inside is multiplied by the Fresnel reflection 
coefficient. 

Following Refs. lHH [HI] the unstable manifold can be 
uniquely computed as a survival probability distribution cal- 
culated from an ensemble of rays starting uniformly in phase 
space having identical intensity. Figure [2] depicts the resulting 
Fresnel weighted unstable manifolds for the lima^on cavity 
using 50 000 rays. Note that (i) in the leaky region, the man- 
ifold is concentrated on very few high-intensity spots. We 
therefore expect a highly directional FFP. (ii) While in the 
case of TE polarization one finds one spot with x > (and 
another symmetry-related one at 5 ^ Smax — 5, X ^ — x)' the 
TM polarization case possesses two of those. 

The unstable manifold in the leaky region directly deter- 
mines the FFP. Mapping the unstable manifold in Fig.[2]to the 
far field by using Snell's and Fresnel's laws (for generaliza- 
tion to curved interfaces, see [33J ) we obtain Fig.O Note that 
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FIG. 1: (color online), (a) Whispering-gallery ray trajectory in the 
circular microdisk, (b) Poincare SOS showing the trajectory (thick 
dots) in phase space; s is the arclength coordinate and x is the angle 
of incidence. Typical trajectories fill a line of constant sin x- The 
critical lines sinXc = ±l/n enclose the leaky region, (c) Chaotic 
ray trajectory in the limafon cavity ([B with s = 0.43. (d) Poincare 
SOS for the trajectory (thick dots) starting above the critical line 
(square) and refractively escaping after only 20 bounces (triangle). 
Small dots show a typical trajectory in the corresponding closed bil- 
liard system. 




0.2 0.4 0.6 O.i 



'.' TE 




- 














... 

- TM 

- 








0.2 


0.4 0.6 0.8 1 







FIG. 2: (color online), (a) Fresnel weighted unstable manifold of the 
limagon cavity for TE polarization. Magnification reveals the differ- 
ences between TE (b) and TM (c) polarization in the leaky region, 
which originate from Fresnel's law. 



the FFP is shown only for the upper half space (0° — 180°), 
the lower half space is given by symmetry. For TE polariza- 
tion, we find directionality around = 0, whereas in the TM 
case additional, smaller peaks occur. The ray in the left up- 
per inset represents one typical trajectory emitting to (j) ^ 
(marked by arrows). The emitting bounce (marked 1, Is is the 
symmetry-related counterpart), the three bounces before and 
the one after (marked 2) are shown. Whereas the trajectories 
are equal for both polarizations, their intensities are different: 
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As visible in the right inset, the rays escaping at 1 hit the Hne 
of the Brewster angle \xb\ = arctanl/n < arcsinl/n. In 
the TE case, transmission is nearly complete and no intensity 
can reach the next bounce 2. This causes the sharp decrease 
in the intensity that is more clearly visible in Fig. Ob). 
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FIG. 3: (color online). Far-field emission pattern for TE (top) and 
TM polarization (bottom) calculated from an ensemble of rays on 
the unstable manifold in Fig. Qb) and (c). Insets illustrate ray dy- 
namics leading to directed emission. The dashed line has strongly 
reduced intensity due to reflection near the Brewster angle. Right in- 
set: 250 rays were started in the rectangular region at 1 and followed 
forwards (region at 2) and backwards (the 3 regions below the critical 
line). Is and 2s mark the symmetry-related bounces to 1 and 2. 

The Eresnel law for TM polarization does not show the 
Brewster angle feature. Therefore, a significant percentage of 
the light is reflected towards bounce 2. Since bounce 2 emits 
into a different direction, an appreciable amount of intensity 
collects in a second far-field peak. 

To summarize up to this point, we have seen that chaotic 
ray dynamics can lead to highly directional emission. How- 
ever, the Q-factor is low since the light rays typically leave the 
cavity very quickly. Does this result of geometric optics carry 
over to the wave dynamics of the electromagnetic field? 

It has been demonstrated that the EFP of optical modes can 
be strongly influenced by the unstable manifold of the under- 
lying ray dynamics lUsl [H, |3ll] . A consequence is that for 
fixed polarization and cavity parameters the EFP is indepen- 
dent on the internal mode structure |30]. This raises the hope 
that the high directionality observed in ray simulations of the 
limagon cavity will survive in wave optics. To this end we 
solve Maxwell equations numerically using the boundary ele- 
ment method [341 . According to the discrete symmetry, even 
and odd modes are distinguished. 

The top panel of Fig. [4] shows near- and far-field pattern of 
a high-Q TE mode. The normalized frequency ft = ujR/c = 
26.0933, c being the speed of light in vacuum, corresponds 
to, e.g., a free-space wavelength of about 900 nm for R = 
3.75 /im. Indeed, as predicted by our ray dynamical analysis 
the mode exhibits directional light emission around = 0. 
The angular divergence of 24° is significantly smaller than 



the values reported for low-Q disks f2T, ^23], and also less 
than in Ref . L241 . Moreover, for fixed polarization and cavity 
parameters, the EEPs of all high-Q TE modes in this cavity 
have similar envelope even though the internal mode structure 
is in general different. This is exemplarily demonstrated in 
the upper panel of Fig.|4]for an odd-parity mode (dashed line) 
which is quasi-degenerate with the even-parity solution (solid 
line). Note that in other types of cavities the even- and the 
corresponding odd-parity solution have in general a different 
EFP 1124]. For the high-Q modes with TM polarization we 
also observe an EFP which is independent of the internal mode 
structure, but a slightly different one, see lower panel of Fig. (H 
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FIG. 4: (color online). Angular dependence of the far-field electric 
field intensity for a TE mode of even (top, solid line) and odd parity 
(top, dashed) with Q = 26.0933, Q ^ 185 000, and a TM mode of 
even (bottom, solid line) and odd parity (bottom, dashed) with Q — 
25.8069, Q ^ 10^; cf. Fig.[3] Insets contain the near-field pattern of 
even-parity modes (middle) and its external structure (right). 

Whereas the ray and wave based EEPs in Figs. [3] and IH re- 
spectively, agree remarkably well, other wave properties seem 
to contradict the ray simulations: (i) the mode does not look 
chaotic but spatially rather well confined, (ii) the cavity Qs 
are too large if compared to the escape rate from the chaotic 
repeller, in fact their values reach, or even exceed, the present 
limit achievable for microdisks with low residual absorption 
and surface roughness | 9]. 

To further investigate the character of these optical modes 
we consider the Husimi projection [ 35], representing the wave 
analogue of the Poincare SOS. From ray- wave correspon- 
dence one would expect that the Husimi projection is dis- 
tributed uniformly over the unstable manifold. However, 
Fig. Oa) demonstrates that the TE mode is localized around 
I sin X I ^ 0.86 and has only exponentially small intensity 
in the leaky region which explains the high Q -factor. A 
closer inspection reveals that the mode intensity is en- 
hanced around an unstable periodic ray trajectory (dots and 
inset in Fig. Oa)) which is part of the chaotic repeller. This 
phenomenon is called scarring [27,1 and has been observed 
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in several kinds of physical systems including microcavi- 
ties iSTllIslllllliO]. Note that the other high-Q modes found 
in this system also exhibit localization along - in general other 
- unstable periodic rays. The Husimi projection of the TM 
mode in Fig. |4] looks similar (not shown). The localization is 
even stronger, leading to a higher Q -factor. 
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FIG. 5: (color online), (a) Husimi projection of TE mode in Fig.|4] 
cf. Poincare SOS in Fig.flld) and unstable manifold in Fig.[3a). The 
lines sinxc = =bl/n enclosing the leaky region are indicated. The 
dots mark the periodic ray trajectory illustrated in the inset. Magni- 
fied Husimi projection in the leaky region for the TE (b) and TM (c) 
mode; cf. Fig.|2tb) and (c) for the ray simulation results. 

Even though the Husimi projection has an exponentially 
small contribution in the leaky region, it is precisely this out- 
going light that determines the FFP. Figures [3b) and (c) show 
the Husimi projection in the leaky region. The convincing 
agreement with the unstable manifold in Figs. Ob) and (c) 
demonstrates its responsibility for the directional emission, 
whereas scarring guarantees the high Q -factor. 



Note that, due to the observed agreement between ray and 
wave simulation, our results are also applicable to larger cav- 
ities. The particular deformation parameter e = 0.43 is the 
optimum value for the localization of the discussed FFPs with 
n = 3.3 but highly localized FFP and high Q-factors can also 
be found for 0.41 < £ < 0.49 (not shown), i.e. fabrication 
tolerances are not crucial. Moreover, we tested that our re- 
sults are robust against variations of the refractive index and 
remain valid for 2.7 < n < 3.9. 

In summary, we have proposed a deformed microdisk as 
a novel cavity design for robust directional light emission 
from high-Q modes. No complicated adjustment of geom- 
etry parameters is necessary, and the emission directionality 
is largely independent from wavelength, cavity size, refrac- 
tive index, and the details of the interior mode structure. The 
latter finding is especially relevant for multi-mode lasing de- 
vices. We trace our, at first sight, counterintuitive results back 
to (i) wave localization along unstable periodic ray trajecto- 
ries ensuring high Q-factors and (ii) escape of rays along the 
unstable manifold of the chaotic repeller leading to directional 
emission. The simplicity of the cavity design allows for easy 
fabrication with a wide range of applications in photonics and 
optoelectronics. The discussed mechanisms are not restricted 
to disk-like geometries but can in principle also be exploited 
for other geometries such as deformed microspheres and mi- 
crotori. 
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